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C^ ' This work investigates a d-p Hubbard model by the n-pole approx- 

I , imation in the hole-doped regime. In particular, the spectral function 

'■^ . A{Lj,k) is analyzed varying the filling, the local Coulomb interaction and 

PlJ ' the d — p hybridization. It should be remarked that the original n-pole 

O . approximation (Phys. Rev. 184 (1969) 451) has been improved in order to 

include adequately the fc-dependence of the important correlation function 
{Sj ■ Si) present in the poles of the Green's functions. It has been verified 
that the topology of the Fermi surface (defined by A(ctj — 0, k)) is deeply 
^ ' affected by the doping, the strength of the Coulomb interaction and also 

'^J . by the hybridization. Particularly, in the underdoped regime, the spectral 

^^ ' function A{uj = 0, k) presents very low intensity close to the anti-nodal 

points (0, ±7r) and (±7r, 0). Such a behavior produces an anomalous Fermi 
surface (pockets) with pseudogaps in the region of the anti-nodal points. 
C*^ ' On the other hand, if the d — p hybridization is enhanced sufficiently, such 

^^ , pseudogaps vanish. It is precisely the correlation function (Sj -Si) present 

in the poles of the Green's functions which plays the important role in 
the underdoped situation. In fact, antiferromagnetic correlations coming 
from {Sj ■ Si) strongly modify the quasi-particle band structure. This 
is the ultimate source of anomalies in the Fermi surface in the present 
approach. 



1 Introduction 

More than two decades after the discovery of the cuprates [1], the theoretical 
description of this phenomenon still represents a challenge for the physicists. 
The fundamental aspects of the cuprates phase diagram can be described in the 
following terms. In the overdoped regime, there is a conventional normal metal, 
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while in the very low underdoped, it is displayed an antiferromagnetic insulator 
(AFI) phase. Between these limits, there is the onset of a superconducting (SC) 
phase. Nevertheless, it is the so-called pseudogap metal [5] above the SC and 
close to the AFI phase which attracts much of the interest. It is widely believed 
that the understanding of such region is the key element to reveal the nature of 
the cuprates. 

The pseudogap phase is characterized by a density of states with low inten- 
sity near the Fermi level leading to an anomalous Fermi surface which is subject 
of intensive investigations. In the experimental field, while the Fermi surface de- 
termined by angle-resolved photoemission spectroscopy (ARPES) is formed by 
disconnected 'Fermi arcs' [3], the quantum-oscillation experiments detect hole- 
pockets enclosing the antinodal points [4]. From a theoretical perspective, some 
general arguments have been proposed to explain the origin of the pseudogap. 
The first one relates it to a preformation of electron pairs [51 IB] . In the second 
one, the pseudogap phase is originated from a hidden broken symmetry [71 [5] . 
However, there is also a proposition which suggests that the presence of short- 
range antiferromagnetic (AF) correlations could be the mechanism responsible 
for the pseudogap [H [10]. For instance, a recent work [11] has indicated that 
antiferromagnetic fluctuations with short correlation length could play a funda- 
mental role in order to understand the apparent disagreement between ARPES 
and quantum-oscillation experiments. Therefore, one can set the question how 
that contributions coming from the antiferromagnetic correlations, could give 
rise to the pseudogap and consequently, the anomalies on the Fermi Surface 
(FS). 

Among the available many body techniques, the n-pole approximation |12[ 
fT3l IM] seems to be quite suitable to such a purpose (see also discussion in 
reference [H]). The structure of the poles of the one particle Green's function 
presents an energy shift (or band shift) in the self energy composed basically 
by three terms, one of them is the correlation function {Si ■ Sj), which can 
produce AF correlations. Actually, recent works have investigated the role of AF 
correlations using the Composite Operator Method (C0M)[TB1[53] which is an 
analytical technique that allows calculating all correlations functions involved 
in the problem, fully self-consistent. As main result, it is obtained that AF 
correlations can be, indeed, a source of a pseudogap. 

It should be stressed that the previous results have been obtained for the one- 
band Hubbard model. One remaining issue is related to the role of hybridization 
on the AF correlations. A recent experimental work [53], using the resonant 
inelastic x-ray scattering (RIXS) technique, has evidenced that the hybridization 
(between the d-orbitals of the cooper Cu and the p-orbitals of the oxygen O) in 
cuprates systems is stronger than that one expected from tight-binding theory. 
However, the one-band models neglect the presence of the oxygen sites. Due 
to the strong correlations at the Cw-sites, the oxygen sites may be occupied 
by holes when the system is doped [H]. Therefore, although the one-band 
models are able to capture many of the most important physical properties of 
cuprates, probably, a model which can incorporate additional details such as 



the hybridization would be a more adequate description of the cuprates. 

In the present work, an extended d — p Hubbard model is treated by a n-pole 
approximation |12| . Therefore, one can investigate how AF correlations could 
affect the properties of the model, mainly focusing on theirs role as a possible 
source of a pseudogap and anomalous FS. It should be remarked that in the 
extended d—p Hubbard model used in the present work the model parameters 
which can affect AF correlations are not only the doping S and the strength 
of local Coulomb interaction U, but also the hybridization. Consequently, it 
would be possible to produce a systematic investigation of how the variation of 
these model parameters as S, U and, particularly, the hybridization, can affect 
the AF correlations and how such correlations influence the features of physical 
quantities such as the spectral function A{w,k). 

In fact, the extended d—p Hubbard model has been previously studied within 
the n-pole approximation |26[ 127] . in which dx2-y2-wave superconductivity has 
been treated following the approach proposed in Ref. [13] ■ In Ref. [17]: it has 
been obtained that the increase of hybridization suppresses the superconductiv- 
ity. This occurs mainly due to the fact that the hybridization has strong effects 
in the band shift which moves the superconducting gap to low energies breaking 
the gap symmetry relative to the zero energy. Furthermore, the hybridization 
broadens and suppresses the density of states becoming unfavored the pair for- 
mation (near the Fermi energy). It is important to mention that in the present 
work, the d — p model includes hoppings to the second-nearest-neighbors which 
improves significantly the uncorrelated band structure considered here. 

It is well known from earlier results jl5) that in the normal state the function 
{Sj -Si) is strongly dependent on the doping and the Coulomb interaction, 
particularly, on the underdoped regime. Moreover, it has also been verified that 
the correlation function {Sj~S~), which behaves like {Sj -Si), is also strongly 
dependent on the d — p hybridization f^P [28 . However, for a given doping 
6, U and tpd produce effects in opposite direction, i.e., while U increases the 
intensity of {Sj ■ Si), f^^ decreases it. The fact that f^^ depletes the correlation 
function {Sj ■ Si) and consequently the antiferromagnetic correlations makes 
the hybridization an adequate physical mechanism to verify whether the AF 
correlations are or are not essential to the onset of the pseudogap. 

It is important to highlight that, in the present scenario, the AF correlations 
are deeply influenced by the momentum structure of the spin-spin correlation 
function {Sj-Si). Indeed, there is a considerable number of works treating about 
this important subject. For instance, in references [TT] [THl [H] the behavior of 
the AF correlations associated to {Sj ■ Si) were studied by using Monte Carlo 
simulations. However, the projection techniques [H] and methods like slave- 
boson [20] and exact diagonalization [21] have also been considered. Actually, 
quite recent results [16] obtained within the composite operator method (COM), 
which is an analytical and fully self-consistent method, have shown that the 
momentum dependence of the charge and spin correlation functions is very 
important to obtain anomalies like the pseudogap in the underdoped regime of 
the hole-doped case. Nevertheless, in the original Roth's procedure [H] and in 



the subsequent works in which the Roth's n-pole approximation is used, it was 
considered tf^ — f^ (where tfj is the hopping between nearest-neighbors) for 
the z nearest-neighbors. Consequently, several correlation functions present in 
the band shift (including {Sj ■ Si)) are constant over all the first Brillouin zone. 
However, such a simplification is unable to capture properly the /c-dependence 
of the band shift Wi , mainly, in the underdoped regime. Thus, in the present 

work, the band shift and, therefore, the (Sj ■ Si) have been evaluated including 
adequately such ^-dependence. 

This paper is structured as follows. Section 2 presents the model and the 
general formalism of the Roth's n-pole approximation 12 . In section 3, the nu- 
merical results are presented and discussed. Finally, section 4 shows a summary 
and some concluding remarks. 



2 Model and General Formulation 

The Hamiltonian model proposed here [211 [27] is given by: 

H = H + H' (1) 

with, 

+ E^^p-^)pIp-+ E ^lAv,.+ E t, [AaVi.+vtdj.) (2) 

where /i is the chemical potential. The term H' , which represents the second- 
nearest-neighbors, is written as: 

^' = E *"4rf..+ E OU- (3) 

The dl„{dia-) and Pi„{Pia) are the creation(annihilation) operators for electrons 
with spin cr in a site i. The quantity U stands for the local Coulomb interaction 
between two d-electrons with opposite spins. The model ^ considers a small 
tfj hopping between d-orbitals and a large tf- hopping between p-orbitals. The 

quantity t^^ stands for hopping between d- and p-orbitals, and A^p — Sp—Sd rep- 
resents the on site-energy difference between the d- and p-orbitals, respectively. 
The parameters tff and t^^ presented in H' represent the hoppings to the second- 
nearest-neighbors for d and p electrons, respectively. The symbols (...) (((...))) 
denote the sum over the first (second)-nearest-neighbors of i. For a rectangular 
two-dimensional lattice, the p-dispersion relation for the first-nearest-neighbors 
is given by: 

e^^ = 2tP[cos{k^a) + cos(A:j;a)] (4) 



where a is the lattice parameter. Considering the homothetic relation [5U1 [SD] 
for the p and d dispersion relations, e'^- = cte^-., where a is a phenomenological 
parameter less than the unity. If the hoppings to the second-nearest-neighbors 



parameter less than the unity, 
are taken into account, 



eZ ^ e'' -. + At^'' cos{kxa) cos{kya) and ei ^ e"^-!^ + At^"^ cos{kxa) cos{kya) . (5) 

In order to obtain the Green's functions within the Roth's method [12], it 
is necessary to define a set of operators {^,1} that describes adequately the 
relevant one-particle excitations of the system. The set of three operators con- 
sidered here is {di„,nf_^dicr,Pia}- These operators must satisfy, within some 
approximations, the relation [An,H]_ = "^^^^ KnmAm, where An are the op- 
erators of the set {An} defined above. This set of three operators results in a 
three-poles approximation for the Green's functions, which in matrix notation 
is written as: 

G(a;) = G(tj)N (6) 

where 

G(w) = N(cjN-E)-^ (7) 

Here, E and N are the energy and the normalization matrices given by 



Enm = ([K,i/]_ , Aj„]^^j) and Nnm = (K,AJ,J(+)) 



(8) 



where [...,...](+)_ denote the (anti)commutator, and (...), the thermal average. 
Considering the set of operators {An} introduced above and the Hamiltonian 
given by equation ([T]), the energy matrix is: 



E 



and the normalization one: 



N = 



ed + ei + Uni^ 


(Erf + ei + U)nt 


ie, + ej + U)ni, 


Unl^+^k-a 


ypd 


nt^VS'' 



dp 



Vr 



nt„Vi^ 



^k 



(9) 



(10) 



rdp(-,rpd\ 



,dpi,pd\ 



with erf = £d — M- The V-^(VT ) are the Fourier transform of t:^{€, ). It 

k k y _ , '-J^ ^ ^J ^ 

is assumed that the system considered here is translationally invariant, then 

^d — „d 



Finally, the quantity Yf:_^ is defined as: 



(11) 



where Wt is the band shift that will be introduced later on. 

One of the most important elements of the Green's function matrix is the 
element 

G^H^ (- - E^^)A~^(-) (12) 

where, 

A^j^iu;) = niM - nt,)[iu-s, - C/(l - nlj - W^J (13) 

and 

X [(w -Sd- ei){u -Sd-U- Wj:J - Unt.iei - M/^rJ]. (14) 

The quantity £'33 is an element of the energy matrix given in equation (|9]). 

Considering the Green's function defined in equation ((T^ , the spectral func- 
tion can be defined as: 

A„{k,uo) = --l^[Gf\u:)]. (15) 



In the real space, the band shift presented in equation (llip is written as 



W,,-. = Wt^_, + W?t„ (16) 



where Wf,_^ is given by equation (51) in Ref. [17]. However, in Ref. p7] . 
the hopping to second-nearest-neighbors has not been considered, therefore the 
band shift Wi presented in equation pT|) must be changed by: 

ni{l-n'i)Wl = h,„+ Y, 4/^-^'^h.2,,+ Y. C,e''-"^h,,, (17) 

which take into account the hopping to second-nearest-neighbors. The quanti- 
ties hia and h2ja are defined as: 

and 

h2j. = B,„ + {Sj ■ So) (19) 

with 

(20) 
The spin-spin correlation function introduced in equation (J19p is given by: 

{Sj -80)^1 {{S+So) + {SiS+)) + (S^S^). (21) 



Particularly, in the paramagnetic state, {S^'S^ ) = (S^ Sq), and the spin-spin 
correlation function {Sj ■ Sq) can be written as: 

{S,-So)^{S+So) + {S^S^) (22) 

where, 



and 



(„d\2 _ O^jcrlT-Oja + Pjam-jcr) 
["'a) 



(24) 



1 - PaP-a 

In Refs. [13l [26l [27l [28], the band shift has been evaluated following the 
original Roth's procedure [12], where ig^ — t^ for the z nearest-neighbors has 

been considered. Consequently, this procedure removes the A:-dependence of the 
correlation functions (mainly {Sj ■ Si) ) present in the band shift. However, the 
momentum dependence of such correlation functions is a fundamental ingredient 
to investigate anomalous properties like the pseudogap in the underdoped regime 
of the Hubbard models [16l[3T]. Therefore, in the present work, the band shift 
defined in equation p7)) is rewritten as: 

«-£ - ;3(T^li:^(«-»W(A (25) 

where the momentum dependence of the correlation functions has been main- 
tained. The e(fc ~ q) is given by 

eik-q)^ ^ t^^e'^-^-^^-^^ + ^ ^u^^(k-q)■R, (26) 

and Fa{q) is given in terms of the Fourier transform of n'^Q^, raja, Cija and jSja 
introduced in equations p^ - (l^^ and defined as: 



<, ^ (4,d,.) - ^E-^-^"'^^'''^^ (27) 

k 

,. = {dlnj^^d,,) = \Y.^-G'^y^^'^'^ (28) 

k 



L 

k 



and 



^Y.^-^^-"^" (30) 



L ^-^ ka 

k 

where J-^^Viuj) = -^ § dujf{Lij)T{ijj), in which /(w) is the Fermi function and 
T{u)) a general Green's function. The Green's functions G and G are obtained 
from the definitions ([6]) and (O , respectively. L represents the number of lattice 

sites. 
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Figure 1: The spectral function A{k,uj = 0) representing the Fermi surface for 
different dopings S. The model parameters considered here are U = 10 eV, 
yj"' = 1.0 eV, t'^ = -0.5 eV, t" = 0.1 eV and keT = 0.01 eV (fcs is the 
Boltzmann constant). 



3 Numerical Results 

In this section, a detailed investigation of the Fermi surface associated with 
the spectral functions for hole-doped regime is done. As a starting point, the 
d— p hybridization has been considered fc-independent [27] (2V^^ )^ = (V-^^V^ ). 
Here, (...) is the average over the first Brillouin zone. The remaining model pa- 
rameters are within reasonable ranges estimated for cuprates j32|. In particular, 
the parameters a == 1, e^ = 0, e^ = 3.6 eV, t^ = —0.7 cV, a = 0.715 and t^P = 
have been kept the same for all results presented here. The homothctic disper- 
sion relation considered in the present work (see equations ^ and ([5])), indeed, 
signify that f^ = at^, therefore f^ ~ —0.5 eV. 

Figure [1] shows the Fermi surface for three different doping levels S, where 
S = 1 — riT (with fiT = n^ + n'L^)- In figure [TJa), 5 = 0.30, and a well defined 
electron-like Fermi surface is observed. In figure[ljb), S = 0.15 and the nature of 
the Fermi surface changes to hole-like. However, it is in the underdoped regime, 
S = 0.07, that the topology of the Fermi surface changes drastically with the 
emergence of a hole-pocket enclosing the nodal point (^, ^). As a consequence, 
due to low spectral intensity, a pseudogap emerges near the antinodal points 
(7r,0) and (0,7r), as shown in figure [Ijc). This result can be better understood 
by analyzing the features of the lower quasiparticle band which is strong affected 
by the antiferromagnetic correlations associated with the correlation function 
{Sj ■ Si). Figure [2^ displays the quasiparticle band for distinct doping levels 6. 
While in the over and moderated doped regimes the quasiparticle bands cross 
the Fermi level near (^, -j) and the antinodal point (0,7r), in the underdoped 
regime the quasiparticle band crosses the Fermi level twice nearer the nodal 
point (f , f ). Such a behavior gives rise to a pocket around (f , f ) (see also figure 
[T};) . On the other hand, as the quasiparticle band does not touch the Fermi level 
near (0,7r), a pseudogap emerges at that region. The kink observed near the 
(tt, tt) point of the quasiparticle band is caused by the strong antiferromagnetic 
correlations associated with {Sj ■ Si), which are maximum in Q= (tt, tt). The Q 
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Figure 2: The quasiparticle bands intercepted by the chemical potential (/x = 0). 
In (a), the bands for U = 10 eV, V^'^ = 1.0 eV, keT = 0.01 eV and three 
different doping levels S. In (b), the bands for 6 = 0.07, V^"^ = 1.0 cV, T = 0.01 
eV and three distinct intensities of Coulomb interaction U. In (c), the bands 
for S = 0.07, C/ = 10 eV, ksT = 0.01 eV and three different hybridizations. The 
figure (d) shows the quasiparticle bands for S = 0.07, U = 10 eV, kq 
and three different temperatures. 



V^^ = 1.0 eV 



U=12.0eV 



(0,0) 




Figure 3: The spectral function A{k,uj = 0) for 6 = 0.07 and VJf = 1.0 eV and 
different intensities of Coufomb interaction. Tfie remaining model parameters 
are identical to figure [TJ 
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Figure 4: The spectral function A{k, a; = 0) representing the Fermi surface for 
different hybridizations Vq = t^'^. The model parameters considered here are 
U = 10 eV, f^ = -0.5 eV, i" = 0.1 eV and ksT = 0.01 eV. 



is the antiferromagnetic wave-vector. 

Figures [2b, Et and [2jl show the lower quasiparticle band for different U, 
Vq and T, respectively. In[2}D and [It, it can be noted that while the Coulomb 
interaction U moves the quasiparticle band near (f ,7r) to lower energies, the 
hybridization V^ moves it to greater energies. Indeed, U and Vq produce 
opposite effects on the pseudogap, i.e., U increases the width of the pseudogap 
while Vq suppresses the pseudogap. Such a behavior of the pseudogap is also 
clear observed in figures [3] and IH 

The features observed in the Fermi surfaces presented in figures [1] |3] and H] 
are understood in terms of the antiferromagnetic correlations associated with 
the spin-spin correlation function {Sj ■ Si) introduced in equations (I19p and 



([22]). Actually, S{k) which is the Fourier transform of {Sj ■ Si), is analyzed 
in the fc-space. Figure [5] shows S{k) at the directions (— 7r,0)-(0, 0)-(7r,0) in 
the first Brillouin zone. In[5][a), with different doping levels S, it is clear that 
S{k) is high doping dependent, mainly in the underdoped region close to the 
antiferromagnetic insulator (AFI) phase. For 6 = 0.30, the chemical potential 
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Figure 5: (a) The Fourier transform {S{ky) of the spin-spin correlation function 
(5,- -5,) for ;/ = 10 eV, Y^^ = 1.0 eV, fc^T = 0.01 eV and different doping levels, 
(b) The S'(fc) for 8 — 0.07 and different intensities of Coulomb interaction. The 
remaining parameters are identical to (a), (c) The behavior of <S'(fc) for 8 = 0.07, 
U = 10 eV, fcsT = 0.01 eV and distinct values oi d — p hybridization, (d) The 
effect of temperature on S{k) for S = 0.07, [/ = 10 eV and V^'' = 1.0 eV. The 



model parameters considered here are t 



-0.5 eV and t 
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0.1 eV. 
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fi is moved to lower energies due to the low occupation nx- As a consequence, 
the Fermi function removes states above the chemical potential leading to a 
suppression of S{k) for k > kp {kp is the Fermi wave vector). The doping 
region, from where the pseudogap emerges, coincides with that in which S{k) is 
very strong. In figure [5jb), the effect of the Coulomb interaction U on S{k) is 
investigated. As it can be verified, U increases the intensity of S{k) mainly at 
the antinodal points region. On the other hand, it has been verified that the d—p 
hybridization V^ ^ acts on S{k) decreasing its intensity (see figure[SJc)). Finally, 
figure[5ljd) shows S{k) for different temperatures. It is interesting to notice that 
the temperature decreases S{k) especially in the region of the antinodal points 
where the pseudogaps are present. This behavior comes from the effect of the 
temperature on the Fermi function, i.e., the Fermi function increases it slope 
near the Fermi level, suppressing S{k) when the temperature is enhanced. If 
the results from figures [5l^a), [Sfb) and[5ljc) are compared with those ones from 
figures [U [3] and m it is clear that the presence of pseudogaps on the Fermi surface 
is directly related to the intensity of S{k) and, therefore, to antiferromagnetic 
correlations. 

It is worth to notice that in the overdoped regime, where the correlations 
associated with S{k) are weakened, the effects of the hybridization are not very 
significant. In the underdoped regime, where mostly sites are single occupied, 
the hopping processes require a double occupied site. As a consequence, the 
electrons tend to stay localized due to the high cost of energy ({/) to double 
occupy a site. Nevertheless, if the hybridization is present, the hopping to a p 
orbital in a site j, via hybridization, can occur. Thus, if the hybridization is 
favored, the occupation of the d-band may decrease depleting the spin-spin cor- 
relations. This scenario allows understanding the effects of the hybridization on 
the antiferromagnetic correlations (associated with {Sj ■ Si)) and, consequently, 
on the Fermi surface topology. 

4 Conclusions 

In the present work, the original two-poles approximation proposed by Roth 
p2] has been improved in order to consider the correct momentum dependence 
of the spin-spin correlation function {Sj ■ Si). The structure of {Sj ■ Si) in the 
momentum space is essential to catch important effects due to antiferromag- 
netic corrrelations which, in the present approach, are the source of anomalies 
as pseudogap and hole-pockets on the Fermi surface. The role of the antifer- 
romagnetic correlations associated with the spin-spin correlation function has 
been investigated in different situations. Initially, it is shown that the hole 
underdoped regime is characterized by the presence of hole-pockets enclosing 
the nodal points and pseudogaps near the antinodal points. The results show 
also that the Coulomb interaction increases the region (in the k space) of the 
pseudogap and decreases the area enclosed by the hole-pockets. On the other 
hand, the Vq ^ hybridization acts in the sense of to decreasing the region (in 
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the k space) where the pseudogap occurs. If the VJ^ ^ is sufficiently high, the 
pseudogap vanishes and an ordinary large Fermi surface is obtained. The sce- 
nario for the Fermi surface described above, can be understood in terms of 
the antiferromagnetic correlations related to the spin-spin correlation function 
{Sj -Si). The hole doping S damps down the intensity of the correlation function 
{Sj -Si). From the overdoped regime up to moderated doping levels {S 0.1), the 
antiferromagnetic correlations are weak and the Fermi surfaces are typical of a 
normal metal without anomalies. In summary, it has been shown that in the 
underdoped regime, where {Sj -Si) is very strong, the pseudogap and the hole- 
pockets emerge. The local Coulomb interaction U increases the {Sj -Si) favoring 
the pseudogap and the hole-pockets. On the other hand, the Vq ^ hybridization 
suppresses the (Sj ■ Si) recovering a large Fermi surface as those ones observed 
in the overdoped regime. Finally, the temperature acts on {Sj ■ Si) decreasing 
it, mainly, at the antinodal points which are the region of the pseudogap in the 
k space. 

To conclude, this work has presented a description for the hole-doped regime 
of the d — p Hubbard model within the n-pole approximation. Particularly, in 
the underdoped regime, it has been shown a route which leads to pockets and 
pseudogap in the Fermi surface. It should be remarked that the correlation 
function {Sj ■ Si) present in the band shift W^ plays an important role in the 
underdoped situation. More precisely, antiferromagnetic correlations coming 
from {Sj ■ Si) strongly modify the quasi-particle band structure. This is the 
ultimate source of anomalies in the Fermi surface in the present approach. 
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